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In quantum dots or molecules with vibrational degrees of freedom the electron-vibron coupling 
renormalizes the electronic charging energy. For sufficiently strong coupling, the renormalized charg- 
ing energy can become negative. Here, we discuss an instability towards adding or removing an 
arbitrary number of electrons when the magnitude of the renormalized charging energy exceeds the 
single-particle level spacing. We show that the instability is regularized by the anharmonic contri- 
bution to the vibron energy. The resulting effective charging energy as a function of the electron 
number has a double-well structure causing a variety of novel features in the Coulomb blockade 
properties. 
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I. INTRODUCTION 

The interplay of mechanical motion and electronic 
properties has been studied widely in molecular physics 
and, recently, in molecular electronics and nanomcchan- 
ical structures. Present-day experimental methods en- 
able the fabrication of nanostructures such as single- 
molecule junctions or quantum dots where a small num- 
ber of electrons couple to one or more vibrational modes. 
Electronic transport in these syste ms exh ibits vibra- 
tional sid ebands at high bias voltages) 1 * 2 * 3 * 4 ^ the shuttle 
instability] 5 * 6 * 7 * 8 * and the Franck- Condo n suppression of 
the low-bias electric conductance! 9 * 10 *^ The shuttle phe- 
nomenon arises from modifications in the spatial con- 
formation driven by the charging events, whereas the 
Franck-Condon effects emerge from the transition rates 
between discrete vibronic states. An additional impor- 
tant consequence of vibrational degrees of freedom is the 
polaron shift which renormalizes the Coulomb repulsion 
of confined electrons and modifies the Coulomb blockade 
physics of the system. At sufficiently strong couplings 
the sign of the effective charging energy U can become 
negative implyi ng an effectively attractive interaction be- 
tween electrons! 12 * 13 * 14 * 15 * 16 * 17 ! The regime of a negative 
effective charging energy is known to be realized in cer- 
tain molecules (an effect known in chemistry as potential 
inversion) ™ It may also be accessible in vibrating quan- 
tum dots when the electron-vibron coupling is sufficiently 
strong 

It has been shown recently that the negative U regime 
opens the possibility for novel features such as an effi- 
cient electron pair tunneling through single molecules^ 
and a charge-Kondo effect P^* An underlying assumption 
of these works is that the negative charging energy does 
not induce an instability of the system. This is in- 
deed the case when the energy gain due to the effec- 
tively attractive charging energy is smaller than the cost 
in single-particle energy due to the finite level spacing, 
when adding (removing) electrons to (from) the dot. This 
paper is devoted to study the opposite regime where 
the magnitude of the negative (renormalized) charging 



energy exceeds the single-particle level spacing. The 
studied system is modeled as a single-electron transis- 
tor (SET) w hose center isl and is coupled to a mechan- 
ical vibration! 9 * 19 * 20 * 21 * 22 * 23 1 It is shown that the system 
becomes unstable towards addition or extraction of elec- 
trons. Assuming that the vibron Hamiltonian contains 
an anharmonic contribution we show that the instability 
is regularized and that the system possesses a well de- 
fined ground state. The effective charging energy as a 
function of the electron number in the dot has a double- 
well structure which, moreover, depends on temperature. 
This is in striking contrast to the usual parabolic charg- 
ing energy and leads to a variety of novel effects. At 
zero gate charge and low temperatures, there is a sym- 
metry between particle and hole-like excitations causing 
the existence of degenerate minima of the effective charg- 
ing energy which for weak anharmonicity are separated 
by a large number of electrons. The average population 
and the electron number fluctuations are highly sensitive 
to the gate voltage which breaks the symmetry between 
the minima. We also show that the transport properties 
of the system exhibit a number of distinctive features 
attributed to the unusual form of the effective charging 
energy. 

This paper is organized as follows. In Sec.[n]we intro- 
duce the model, demonstrate the origin of the negative-?/ 
instability, and state the precise conditions under which 
it takes place. In Sec. |III| we explore the equilibrium 
properties of the instability, including the effects of the 
vibron nonlinearity. Nonequilibrium (transport) proper- 
ties in the presence of a bias voltage are discussed in Sec. 
UV1 We conclude in Sec. El 



II. MODEL AND NEGATIVE U INSTABILITY 

We are considering a quantum dot or a single- molecule 
junction coupled to a single vibrational degree of free- 
dom. The dot is modeled as a SET where the number 
of electrons on the island couples to a vibron mode, see 
Fig. [T] The system could be realized, for example, by 
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FIG. 1: Equivalent circuit of the quantum dot. Capacitances 
Ci and C2 model high resistance tunnel junctions which al- 
low the island to exchange electrons with the reservoirs, C g 
represents a gate capacitance. The central part is coupled to 
a vibrational mode. 



case a Ml <C 1, is crucial for the stability of the system 
for large couplings A. To demonstrate this let us first 
assume a = and examine the phonon-induced modifi- 
cation of the charging energy. Since the electron number 
couples linearly to the phonon position, it is possible to 
introduce a shifted phonon coordinate by completing the 
square in the phonon Hamiltonian. As a result one elim- 
inates the coupling term and generates an extra term 
AU(N) quadratic in N which can be thought of as a 
modification of the charging energy 



U cS (N) = U(N) + AU(N) = 



(Ne 



2C 



la)2 -N 2 x 2 m. 

(5) 



It is important to recognize that the well-known polaron 
shift AU(N) gives rise to a negative contribution. The 
overall prefactor of the N 2 -term will be negative if 



a suspended center clement between large metallic leads 
enabling tunneling between the dot and the leads. The 
electronic part of the Hamiltonian for a SET is 



L.R.C 



-j c j c j 



H T + U(N), 



(1) 



where the first term represents the single-particle ener- 
gies of the leads and the center part, the second term de- 
scribes tunneling between the leads and the center, and 
the last term is the relevant charging energy of the sys- 
tem. The tunneling Hamiltonian and the charging energy 
take the standard forms 



and 



U(N) 



(Ne - Q g ) 2 
2C 



(2) 



(3) 



respectively. Here, N = X)c c j c i — Nq is the number 
of extra electrons on the island, C = C\ + C2 + C g , and 
Q g = ViCi+V 2 C2+V g Cg is the total induced gate charge. 
The number No corresponds to the positively charged 
lattice ions on the central island. The full Hamiltonian 
of the system H = H e + iJ p h also includes the vibrational 
contribution 



H 



P h 



2M 



-Mfl 2 x 2 + a(2Mn 2 ) 2 x i + \hflN- J ^- (4) 
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where we have included the electron-vibron coupling 
XMlNx/l osc with Z osc = \J H/2MQ, into the phonon 
Hamiltonian. The Hamiltonian Q also includes an an- 
harmonic x 4 term whose strength is controlled by the 
parameter a. The anharmonic term becomes compara- 
ble to harmonic effects when the phonon displacement is 
of the order of (aMl)~ 2 times the oscillator length Z osc . 
The inclusion of this term, even in the weakly nonlinear 
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2C 



X 2 Ml < 0, 



(6) 



meaning that the effective potential favors large absolute 
values of N. Defining the quantity Ec — e 2 /2C, the 
negative charging energy condition ^ takes the form 



A > 



(7) 



Consequently, in metallic quantum dots with a vanish- 
ingly small single-particle level separation, the energy can 
always be lowered by adding or removing an arbitrary 
number of particles when condition ^ holds. Thus the 
system does not possess a ground state and is unstable. 
In dots whose orbital energies due to confinement are 
not negligible compared to the charging and vibron ener- 
gies, the single-particle energy cost should also be taken 
into account. To analyze this, we introduce equidistant 
single-electron orbit als 



nir 

E n = hv F — = An, 
I'd 



(8) 



where vp is the Fermi velocity of the dot, Ld is the effec- 
tive size of the electronic confinement and A = hvpir/Ld 
is the level separation. Expression ^ is appropriate in 
the vicinity of the Fermi energy and holds exactly for 
a one-dimensional Dirac spectrum. The energy cost of 
adding or removing TV electrons to or from the lowest 
available orbitals (ignoring the spin degeneracy) is 



N 



N(N+ 1) 



(9) 



n=l 



Adding the single-particle contribution to the effective 
potential ^ we get 



U eS (N) = A 



N(N + 1) 



+ E C (N - Qg/e) 2 - N 2 \ 2 Ml 



A 



+ E C )(N- Qg/ef - N 2 \ 2 Kl + const, (10) 
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where the second line follows by shifting the gate charge 
Qg — > Qg + 2^(1 + Qg)- The constant in Eq. (10) is 
independent of N and can be dropped. Thus the func- 
tional form of the effective potential remains invariant, 
and the finite level separation just renormalizes charging 
energy and gate charge. The system is still fundamen- 
tally unstable provided that 



A > 



E c + A/2 



(ii) 



In the spin-degenerate case A should be replaced by A/2. 
Alternatively, we can write this condition as Ec< — X 2 hfl < 
—A/2, i.e., we find an instability when the magnitude of 
the renormalized charging energy becomes larger than 
the level spacing. 

Result (111 shows that a constant (or decreasing) level 



spacing is not adequate to stabilize the system at suffi- 
ciently large couplings. The instability tends to change 
the electron number by increasing the displacement of 
the phonon mode. Therefore it is natural to assume 
that at some point when the displacement becomes large 
the harmonic approximation for the vibron mode breaks 
down and anharmonic effects become significant. It turns 
out that a generalization of the above analysis to include 
a nonzero anharmonic term (a > 0) in Eq. Q always 
produces a well-defined ground state. The effective po- 
tential for electrons becomes temperature dependent and 
exhibits a variety of novel features. 

Our discussion has a close relation to familiar stabil- 
ity conditions in Fermi-liquid theoryP^In fact, it is well 
known that the Fermi liquid remains stable even when 
the Landau interaction parameter F becomes negative 
as long as Fq > —1. Since Fq measures the interaction 
strength in units of the density of states, this is a pre- 
cise analog of the stability condition (111 expressed as 
(E c - A 2 M7)/(A/2) > -lP 



III. THERMAL EQUILIBRIUM PROPERTIES 

A. Effective potential 

In this section we study effects of the anharmonic 
phonon term on the negative-C/ instability, focusing on 
thermal equilibrium properties of the junction. Diago- 
nalizing Q one obtains phonon eigenvalues {Ej(N)} as 
a function of the electron number N and the island par- 
tition function becomes 33 



Z = £e X p{-/3[f/(iV) + ^(iV)]} 



N,j 



^exp{-/?C/(7V)}^exp{-/3i?,(7V)}. (12) 



N 



In Eq. ( |12| we have assumed that the electron tunnel- 
ing is weak and its contribution to the partition func- 
tion is negligible compared to the charging energy.^ 



Defining a phonon partition function as Z p h(N) — 
J2j ex P the partition function takes the 

form 



Z = Z ph {N = 0) 
x^exp{-/3[[/(A0 



N 



ft 



In (Z ph (N)/Z ph (N = 0))] 



(13) 



From expression (13 1 one can identify the effective elec- 
tronic potential 

U cS (N, T) = U{N) - hn [Z ph (N)/Z ph (N = 0)] . (14) 

In the absence of the electron-phonon interaction, Eq. 



( 14 1 reduces to the ordinary charging energy, and in the 
absence of the anharmonic term (a — 0) it is tempera- 
ture independent and coincides with expression (J5J). In 
the anharmonic case the polaron shift depends on the 
phonon state and thus the effective potential depends on 
the phonon temperature. In general, an evaluation of 
Eq. ( 14 1 requires one to diagonalize ^ numerically and 
to compute the phonon partition function. The behavior 
of the system is then determined by the relative strength 



of the energy scales Ec- hCl, 



and the magnitude of 



dimcnsionlcss electron-phonon coupling A. 



B. Analytical considerations 

Analytical expressions for the effective potential can 
be obtained for weakly unstable systems at zero temper- 
ature when considering the vibronic degree of freedom as 
classical. As this already illustrates some of the essential 
physics, we first consider this case before presenting more 
general numerical results for the fully quantum mechan- 
ical model. 



Starting with Eq. (14 1, we can write the effective po- 



tential at zero temperature as 



U cS (N) = U(N) + E (N) - Eq(0), 



(15) 



where Eq(N) denotes the ground state energy of the 
phonon Hamiltonian H p h for excess charge N. In the 
classical limit, the ground state energy is simply given 
by the minimum of the potential energy. Noting that 
-Eo(0) vanishes, we obtain 

U cS (N) = U(N) - \ 2 nnN 2 

+ mm{-Mn 2 (x + 2XN£ osc ) 2 + a(2Mfl 2 ) 2 x 4 }. 
x 2 

(16) 

For a weakly unstable system near zero gate charge, we 
expect that the effective potential has a minimum for 
small (but nonzero) N. Thus, we can neglect the quartic 
term when determining the position of the minimum in 
Eq. pi 



x 



-2MIN. 



(17) 
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It is straightforward to include corrections to this expres- 
sion and find that this approximation is valid as long as 
X 2 \Ec — \ 2 Ml\/Ml <C 1, i.e., as long as the instability is 
only weakly developed. With this expression for xq, we 
find for the effective potential 



C/ eff (iV) = U(N) - X Z KIN' Z + 16a(hn)' A N 4 



(18) 



Focusing on the case of zero gate charge, we can now 
compute the number of electrons Nq which are enter- 
ing the quantum dot as a consequence of the instability. 
Minimizing the effective potential in Eq. ( 18 1 , we find 



\E C - X 2 f£l\ 
32a(hn) 2 



(19) 



Note that there are indeed degenerate minima at N = 
±Nq for zero gate charge. 

Once the instability develops, we can define an effec- 
tive charging energy which decribes the curvature 
of the effective potential in the vicinity of the minima. 
Expanding Eq. (18 1 around the minima, we find 



x 2 m\ 



(20) 



This sign reversal of the effective charging energy induced 
by the instability has important consequences. While a 
negative renormalized charging energy leads to pair tun- 
neling in the absence of the instability, 9 the positive effec- 
tive charging energy E^? implies that transport will be 
dominated again by single electron sequential tunneling 
processes once the instability develops. It is also inter- 
esting to note that the Coulomb blockade is significantly 
weakened near the instability where many charge states 
have similar energies. 

We close this section with a remark on the validity 
of the classical approximation for the vibrons employed 
in this section. As can be readily seen from Eq. (15 1, 



the quantum-mechanical zero-point energy of the vibron 
mode cancels out from this expression as long as the vi- 
bron remains harmonic, implying that the zero-point en- 
ergy is independent of the electron number N. Thus, 
quantum corrections to our classical discussion are pro- 
portional to the anharmonicity and our classical expres- 
sions are an excellent approximation to the results of a 
fully quantum mechanical calculatio n as l ong as the an- 
harmonicity remains small (see Sec. Ill C for a compari- 
son with numerical results). 



C. Numerical results 

A solution of the full quantum problem for arbitrary 
parameter values requires a numerical approach. We now 
focus on the results of such a calculation. Figure [2] illus- 
trates the dependence of the effective potential on the 
electron-phonon coupling strength and the charging en- 
ergy. As the coupling strength is increased the system 
exhibits a crossover between the 7 > regime and the 




FIG. 2: Effective potential as a function of the electron 
number. Left: the coupling-dependence of the potential, 
T = 0.1 Ml, a = 0.01, E c = fifi, Q g = and from bottom to 
top A 2 = 3, 2.5, 2, 1.5, 1. Right: the effective potential cor- 
responding to different charging energies, a — O.Ol(ftO) -1 , 
A 2 = 2, T — O.lhQ, Q g = 0, and from bottom to top 
E c = 0.5, 0.8, 1, 1.5, 2 x HQ. 



negative-?/ regime 7 < 0, signalled by the formation of 
a double- well structure. Decreasing the (bare Coulomb) 
charging energy results in qualitatively similar behavior 
as increasing the coupling strength. 



u eS /nn 






J 

N 



10 -5 




FIG. 3: Effective potential as a function of the electron 
number. Left: Curves illustrate the temperature depen- 
dence of the potential corresponding to parameters Ec = tiO., 



a = 0.01(ftfi)~ 



A = 2, 



= at T = 0.1, 3, 5, 10 x HQ 



(from bottom to top ). Right: The same as the left figure but 
at finite gate charge Q g — e/2. 



Figure [3] shows the temperature dependence of a typi- 
cal effective charging curve U eS (N,T) in the negative-C/ 
regime 7 < 0. At low temperatures the potential ex- 
hibits a double-well structure whose details depend on 
the particular parameter values. The shape of the po- 
tential can also change qualitatively as the temperature 
is increased: the double well structure present in Fig.[3]is 
eventually deformed to a single-well potential. The tem- 
perature dependence of the potential is a consequence 
of the nonlinear phonon interaction. In the asymmet- 
ric case (i.e., at nonzero gate charge Q g ^ 0) there ex- 
ists a metastable energy minimum whose decay towards 
the global minimum requires multiple electron tunneling 
(roughly 10 electrons for the parameter values). Figure 
[1] (left) shows a comparison between the full quantum 
mechanical solution obtained numerically and the clas- 



sical result obtained in Se. IIIB As argued above, we 
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FIG. 4: Left: Comparison of the classical and quantum me- 
chanical effective potentials. Solid curves correspond to the 
quantum mechanical calculation with parameters Ec = hfl, 
a = O.Ol(fifi) -1 , Q g = 0, T = and A 2 = 2.5 (blue), A 2 = 2.0 
(green), A 2 = 1.5 (yellow). The dotted lines correspond to the 
classical approximation obtained by minimizing the phonon 
potential energy. Right: Electron number corresponding to 
the minimum of the right potential well as a function of the 
nonlinearity, E c = Hil, A 2 = 2, Q g = and T = 0.1M1 
(solid), T = 3hQ (dashed), T = 5hQ (dotted) and T = WhQ 
(dash-dotted). 



find that the classical approximation is quite accurate 
for small anharmonicity. The wells are closer at stronger 
anharmonicity and higher temperature, as indicated by 
Fig. |4] (right). 




FIG. 5: Average dot population as a function of the gate 
charge. The left figure corresponds to parameters T = 0.1 HQ, 
a — O.Ql(ftfl) -1 , Ec = Ml and from bottom to top A 2 = 
0, 0.2, 0.8, 1, 1.2 and 1.5. At vanishing coupling (A = 0) we re- 
cover the usual Coulomb steps, in the strong coupling we have 
a rapid increase at low gate voltages. The right figure clarifies 
the small gate charge region for couplings A 2 = 0.8, 1, 1.2, 1.5 
(from bottom to top). 



The average number of electrons as function of the 
gate charge is presented in Fig. [5] In the absence of the 
electron-phonon coupling and at low temperatures, the 
average number (N) exhibits the usual Coulomb staircase 
behavior. As the coupling is increased, this is gradually 
transformed into a new dependence reflecting the double- 
well nature of the effective potential. The double-well 
structure leads to a rapid increase of the population at 
low gate charges. The reason for this is that the gate 
acts as a symmetry breaking field and the values Q g /e > 
k^T/Ec are sufficient for the system to prefer one of the 
two nearly degenerate wells. 



This physics is also clearly reflected in the average 
number (N) as function of temperature, as shown for 
various gate charges in Fig. [6] For finite positive gate 
charge the zero temperature average value corresponds 
to the minimum of the right well. As the gate charge in- 
creases, the position of the minimum shifts and electrons 
are added to the dot in discrete jumps of one electron. 
At higher temperatures, the system is no longer trapped 
in the right well so that the average number of electrons 
decreases rapidly. However, the average number satu- 
rates at a finite value at high temperatures, reflecting 
the gate-induced asymmetry of the potential. 
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FIG. 6: Average dot population as a function of temperature. 
The left figure corresponds to parameters A 2 = 2, a = 0.01, 
Ec = 1 and from bottom to top Q g /e = 0.05, 0.1, 0.2 0.4 and 
0.5. The right figure clarifies the low-temperature region. 




FIG. 7: Root mean square fluctuation as a function of temper- 
ature. The left figure corresponds to parameters Ec = ftfi, 
a — 0.01(ftf2) _1 , A 2 = 2 and from bottom to top Q g /e = 
0.5, 0.1, 0.01, 0. At vanishing gate the low-temperature fluc- 
tuations saturate at finite value, otherwise they will go to 
zero. The right figure clarifies the low-temperature region. 

Charge fluctuations also exhibit interesting behavior. 
In the absence of the electron-phonon interaction the 
electronic potential is U(N) — E C (N — Q g /e) 2 , so the 
mean square number fluctuation 

AiV 2 = (iV 2 ) - (AO 2 (21) 

is proportional to temperature, AA^ 2 oc T. The situation 
is very different when the electrons couple to an anhar- 
monic vibron mode. First of all, the zero-temperature 
number fluctuations do not vanish at zero gate charge, 
Q g = 0, but saturate to a finite value, see Fig. [7] This 
happens because there exist symmetric minima on both 
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sides of the origin. The curve also exhibits a weak dip 
at finite temperatures since charge fluctuates more fre- 
quently towards values smaller than the charge at the 
minima. Eventually fluctuations grow due to the popula- 
tion of higher energy states. At finite gate values Q g ^ 
the fluctuations are drastically modified at low tempera- 
tures. The symmetry between the minima is broken and 
only states close to the preferred minimum are populated. 
When there is a single global minimum, the fluctuations 
vanish at zero temperature. At finite gate charge the 
curves approach to the zero gate curve at finite temper- 
atures. At small gate charges, even a weak thermal exci- 
tation is sufficient to restore an approximate symmetry 
between the minima, leading to a sharp increase of the 
fluctuations from zero to the symmetric value at Q„ = 0. 



IV. TRANSPORT PROPERTIES 



suppressed tunneling current Iq 



Ml/eR t . At tem- 



peratures T > Ec /2 current is a slowly varying function 
of the gate charge. Below this temperature it gradu- 
ally starts to show signatures of Coulomb-like oscilla- 
tions which become pronounced at low temperatures. Al- 
though this is in qualitative agreement with the standard 
Coulomb blockade results in absence of the vibron in- 
teraction there are also some essential differences. All 
the curves exhibit a slowly decreasing tendency and the 
amplitude of current oscillations grow as a function of 
gate charge. The strict periodicity of the usual Coulomb- 
blockade conductance is broken by the nonlinear phonon 
interaction which does not leave the effective charging 
spectrum invariant as Q„ is increased by multiples of the 
electron charge. Figure M (right) shows also that, within 
a reasonable accuracy, the current is a linear function of 
the applied voltage at low bias. 



So far we have explored consequences of the negative- 
U and the anharmonic phonon effects on equilibrium 
properties. In this section we focus on nonequilibrium 
characteristics of the system in the negative- U regime. 
To simplify the analysis we assume a continuous den- 
sity of states of the dot and that the energy scales sepa- 
rate Ml » Ec so that only the lowest vibrational state 
\Eo(N)) is relevant at low temperatures and bias volt- 
ages. As implied by considerations in Section III and 



the estimate for the effective charging energy Eq. ( 20 1 



close to the instability the curvature of the potential is 
weak and one expects the sequential tunneling processes 
to dominate the transport. In the lowest order of the tun- 
neling coupling transport properties are described by a 
rate equation with Golden- Rule transition ratesP^ Rates 
for the processes in which the initial dot state contains 
N electrons and where the electron number changes by 
one due to the tunneling through the left/right junction 
are 



N = 7 L/R / (U cS (N ± 1) - U cS (N) ± ML/R ) , 

(22) 

|(g (jv±i)|i; (jv))| 2 , 

e 2 R t 



r 



where fix 



eP'-X 
47re 2 |t| 2 ^ L / R i/ c 
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L/R 



1/Rt = — ^ L/R ^ c ■ In the positive-C/ regime in the 
absence of the anharmonic term the matrix element 
\{E (N ± 1)\E Q (N))\ 2 leads to the usual Franck-Condon 
suppression of the conductance given by \(E (N ± 
1)\E (N))\ 2 = e- x \ The probability P(N) of having N 
extra electrons on the dot follows from the detailed bal- 
ance condition r n-i,nP(N) — T NiN _iP(N — 1), where 
the rates T are given by the sum of the corresponding 
left and right lead rates. The stationary current through 
the SET can then be calculated from the expression 



I = 



N 



r 



N+1,N 



-T 



N-1,N 



P(N). 



(23) 



In Fig. [8] we have plotted the current as a function of 
the gate charge at different temperatures and bias volt- 
ages. The current is normalized to the Franck-Condon 




FIG. 8: Current as a function of the gate charge. Left: Differ- 
ent curves correspond to temperatures (from top to bottom) 
T = 0.2, 0.15, 0.1, 0.07, 0.05, 0.03, 0.02 x htt, other parame- 
ters being E c = 0.2 M7, a = 0.05(ftn) _1 , A 2 = 0.5, and fi L = 
— fiR — O.lhQ/e. Right: Curves correspond to voltages (from 
top to bottom) yL L = -fi R = 0.1,0.08,0.06,0.04,0.02 x ftfi/e, 
other parameters being Ec = 0.2 hQ, A 2 = 0.5, a — 0.05, and 
T = 0.05 hfl. 



P ( N ) h 



N 




FIG. 9: Left: Probability distribution of the dot population 
corresponding to parameters Ec = 0.2HQ, a = 0.05 (FiQ) -1 , 
A 2 = 0.5 and fi L = -fx R = O.lhQ/e, Q g = at T = O.lffi. 
Right: The effective potential corresponding to the same pa- 
rameters. 



It is noteworthy that the current is not completely 
blocked even at relatively low temperatures T ~ 0.1 Ec 
where the usual Coulomb-blockade current is completely 
suppressed around integer gate charges. The reduction 
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of the blockade can be qualitatively understood by con- 
sidering the estimate ( 20 1 , which indicates that in the 



vicinity of the instability the value of the effective charg- 
ing is reduced. With parameters corresponding to Figs. [8] 
and [9] the minima of the double- well potential are flat- 
ter than the minima corresponding to the bare electronic 
potential in the absence of vibrations and, as illustrated 
in Fig. [9j there exists two (and at weaker nonlinearity 
more) almost degenerate states in each well. Transitions 
between these closeby states at the bottom of the wells 
enable current to flow below temperatures where the bare 
electronic energy cost would block it. The estimate for 
the effective charging ( 20 1 is valid only in the vicinity of 
the instability and by increasing \Ec — XMl\ one even- 
tually enters deep in the negative- U region where the 
curvature at the potential minimum exceeds Eq- In this 
region the blockade is stronger than the usual Coulomb 
blockade and one has to consider higher-order tunneling 
processes. 

In Fig. 10 we have plotted current vs. gate charge 
with a variable phonon nonlinearity. Small changes in 
the strength of the nonlinearity change energy differences 
of the nearby states at the potential minima leading to 
phase shifts in the current oscillations. Since the current 
is an even function of the gate charge it exhibits a cusp 
close to Q g = depending on the phase of the current os- 
cillations. The most pronounced effect is achieved when 
the current oscillations jump by half a period at the ori- 
gin. At finite temperatures the cusp is always smooth and 
becoming sharper when temperature decreases. The rea- 
son for the existence of the cusp is the symmetry breaking 
between the two minima of the effective potential at fi- 
nite Q g . Even small gate values localize the system in 
one minimum leading to observable effects in the current 
through the structure. 
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FIG. 10: Current as a function of the gate charge at differ- 
ent phonon nonlinearity strengths. Left: The different curves 
illustrate cases a = 0.058 (fifi)" 1 ( green), a = 0.05 (hfl) 
(red), a = 0.045 (fin) -1 (blue) and a = 0.04 (M7)- 1 (black), 
other parameters being Ec = 0.2fif2, A = 0.5, T = 0.05 HQ. and 
= — fin — 0.05liQ/e. Right: Same quantities in the vicin- 
ity of Q g — 0. The curves correspond to a = 0.056 (hQ)^ 1 
(black), a = 0.053 (M7) -1 (green), a = 0.05 (hQ.)' 1 (red) and 
a = 0.047 (ftf2) _1 (blue), other parameters as in the left plot. 



V. CONCLUSIONS 

In this paper, we studied effects of strong electron- 
phonon interaction in quantum dots with a vibrational 
degree of freedom. At sufficiently strong couplings the 
vibron-induced polaron shift overcomes the charging en- 
ergy cost and the effective potential for electrons fa- 
vors large charging. The instability towards an arbi- 
trarily large electron population on the dot is regular- 
ized by the anharmonic contribution to the phonon en- 
ergy. The effective potential differs qualitatively from the 
usual Coulomb repulsion, leading to characteristic modi- 
fications of the low temperature Coulomb blockade prop- 
erties. Signatures include a rapid change in the average 
and the fluctuations of the electron number as function 
of the gate charge in the neighborhood of Q g = as 
well as the temperature dependence of these quantities. 
Moreover, transport properties also show a number of 
unique signatures that can be used to characterize the 
negative-?/ instability and the nonlinear phonon effects. 

It is interesting to check whether in addition to 
molecules, an effectively negative U could also be 
achieved in nanoelectromechanical systems. As shown 
by Eq. (Ill, this requires an electron- vibron coupling (as 



measured by Xhuj) which is of the same order of mag- 
nitude as the charging and the single-particle energies. 
It turns out that this condition rules out currently ex- 
isting suspended electron beam lithography samples. At 
the same time, the situation is more favorable (though 
still marginal) for suspended carbon-nanotube devices, 
for which strong electron-vibron coupli ng ha s been ob- 
served in several recent experiments! 26 ' 27 ^ The most 
likely vibron mode to cause a negative U in this system 
is the radial breathing mode with frequency u) ~ c/Lj_- 
(Here, c denotes the velocity of propagation of acous- 
tic phonons in graphene and L± is the circumference 
of the nanotube.) The charging energy and the single- 
particle level spacing are of the same order (assuming 
a "fine-structure constant" e 2 /Hvf ~ 1 as appropriate 
for graphene) so that in order to reach the negative U 
instability, the electron-vibron coupling A must exceed 
[{vf / c){Lj_/ L)] 1 / 2 for a nanotube of length L. Using 
vf/c ~ 10 2 , L± ~ lnm, and L ~ 1/im, we find that the 
critical coupling strength is A c ~ 0.1 — 1. Unfortunately, 
this coupling constant is relatively poorly understood at 
present. However, both experiment^ 6 -' and theoretical es- 
timates for one possible coupling mechanism 29 suggest 
that it is of the same order as A c p2 While these estimates 
do not allow us to claim that the instability discussed in 
this paper could be observed in carbon-nanotube devices, 
they do suggest that realizing a negative U in nanoelec- 
tromechanical systems is certainly conceivable. 
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Tunneling is still important in the sense that in the equi- 
libration process the island can change its charging state 
by exchanging electrons with the leads. 
Strong electron-vibron coupling has been observed re- 
cently, e.g., in suspended carbon nanotube quantum 

dots flnaSkan 

Quantity (Ec — A 2 ftf2)/(A/2) corresponds to Fo apart the 
difference of factor two which can be readily traced to triv- 
ial differences in definitions. 

For simplicity, we present our analysis in the case A = 0. 
A nonzero A can be readily included by a simple renor- 
malization of Ec and Q g . 

According to Ref. 1291 A for the breathing mode scales as 
L~ 3 which means that decreasing the critical coupling 
strength by increasing the length of the tube does not to 
improve the possibility of reaching the negative-f/ regime. 



